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Abstract 



We study the effect of the envhonment on the process of the measurement 
of a state of a microscopic spin half system. The measuring apparatus is a 
heavy particle, whose center of mass coordinates can be considered at the 
end of the measurement as approximately classical, and thus can be used as 
a pointer. The state of the pointer, which is the result of its interaction with 
the spin, is transformed into a mixed state by the coupling of the pointer to 
the environment. The environment is considered to be a gas reservoir, whose 
particles interact with the pointer. This results in a Fokker-Planck equation 
for the reduced density matrix of the pointer. The solution of the equation 
shows that the quantum coherences, which are characteristic to the entangled 
state between the probabilities to find the pointer in one of two positions, 
decays exponentially fast in time. We calculate the exponential decay function 
of this decoherence effect, and express it in terms of the parameters of the 
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I. INTRODUCTION 



In recent years there has been a considerable interest in the effect of the environment on 
the process of quantum, measurement. The measurement of physical variable of a microscopic 
system - a system which is described by quantum theory - is usually realized by coupling the 
system to a "classical meter". The meter itself does not have to be macroscopic, however 
some of its physical variables, say the position of its center of mass, can be considered as 
approximately classical. If these variables are chosen to be correlated with those of the 
system, their readings are related to the measurement of the microscopic properties of the 
system. When the measurement process is efficient there is one to one correspondence 
between the macroscopic readings of the meter and the microscopic quantities, which are 
set to be measured. A particle, whose center of mass position serves as a meter, will be 
identified as a pointer. A theory that starts from a model Hamiltonian, which includes 
only the system, the pointer, and the coupling between them, cannot account for the real 
measurement process. The reason is that the classical meter is significantly infiuenced, on 
the time scale of its motion, by the environment around it. Even though the dynamics of 
the combined entity - the system, the meter, and the environment - can be described by 
quantum mechanics, the environment (bath) plays a distinctive role in the process. 

The role of the environment is formally expressed in the quantum description of the 
experiment in that the system ends up in a mixed state, rather then in a pure state, which 
is the outcome of a "pure" unitary transformation. While a unitary transformation leaves 
the entire entity, system, meter and bath, in an entangled state, so that the pointer's final 
positions are not definite, but rather display interference, the environment generates deco- 
herence, and induces definite final readings of the meter. The fundamental nature of the 
quantum theory makes the predictions of the readings probabilistic, while the decoherence, 
which is caused by the nature of the environment, is responsible for the definite positions of 
the pointer. 

That decoherence is induced by the environment was recognized long ago, as was pointed 
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out explicitly almost twenty years ago by Zurek . This point was previously noticed by 
Zeh [0 and Kiibler and Zeh |^ in the early seventies. Soon later Caldeira and Leggett , 
and Walls, Milburn and Collet [p,^ have shown, using a harmonic oscillator as their micro- 
scopic quantum system, to be measured, and an ensemble of many harmonic oscillators to 
model the environment, that indeed the environment introduces decoherence to the Density 
Matrix of the quantum system. In Ref. an initial pure state of two spatially separated 
Gaussian wave-packets of the microscopic system - a harmonic oscillator - are turned by the 
bath, in an extremely short time, into a well defined mixed state of the two wave-packets. 
Ref. starts from a pure state with a superposition of two separated pure coherent states 
of a harmonic oscillator, and ends up, again due to the interaction with the bath, in a mixed 
state of the two coherent states. The quantum interferences between the two wave-packets, 
or the two coherent states, fade away in an extremely short time. Kiibler and Zeh H] and 



later Zurek et al. [T^l have shown that the coherent states of a harmonic oscillator in a bath 
are somewhat unique in their decoherence, being almost classical. Collett ||Tl| has extended 
the investigation of decoherence using density-matrix calculations for some simple open sys- 
tems. Following the ingenious experiments of the Ecole Normal Superieure Group Paz 
and Zurek |T3] have investigated the decoherence in the limit of weak interaction between 



the "meter" and the bath. Recently the decoherence of a harmonic oscillator, playing the 
role of a meter measuring a spin half system, has been studied by Venugopalan [jl4|. Their 
model for the bath was the standard one 0, namely an ensemble of harmonic oscillators. 

In the present paper we revisit the decoherence problem, and study it more in the context 
the "measurement" theory in a somewhat different setup. We consider a simple model to 
describe a combined system, which is made of three coupled parts, a microscopic atomic 
system, a measuring apparatus, and a bath. While the standard model for the bath is an 
ensemble of independent harmonic oscillators, we consider here a gas reservoir, namely an 
ensemble of independent particles, which interact with the meter. The initial state of a spin- 
half atom is set to be "measured" by the final position of a massive particle, which is, in 
turn, in contact with the reservoir. This combined system is supposed to simulate the entire 



process of a "real" measurement. An environment of a gas seems to be more appropriate for 
the description of the motion of a particle-pointer. Furthermore, by changing the density 
of the gas one can control the time of decoherence. The atom interacts for a very short 
time with the meter, and the state of the two systems is becoming entangled. The two 
macroscopically separated positions of the pointer, which are now quantum-mechanically 
correlated with the initial state of the spin, interfere with each other. The interaction of 
the meter with the particles of the bath introduces decoherence, namely is destroying the 
interference, and making the pointer's readings distinguishable. The development in time 
of the entire system is described by a unitary transformation, which originates probability 
into the dynamics by quantum mechanics. This evolution in time is conveniently described 
by the density-matrix of the entire system. In the treatment of the measurement process, 
a pure state is transformed after integrating out the environment into a mixed state, and 
the notion of the collapse of the wave function is never being introduced. The probabilistic 
nature of microscopic dynamics, which is built into the quantum description, is expressed, 
at the end of the process, in terms of the two potentially definite positions of the pointer. 

The crucial role of the environment is clearly expressed in the stage where the density- 
matrix of the combined system of the pointer and the bath is traced over the degrees of 
freedom of the gas particles of the reservoir. Since the exact initial conditions of the bath 
cannot be completely prescribed, the unitary transformation, which controls the motion in 
time, does not account for the complete description of the entire process. The point is that 
even if we could theoretically prescribe the state of the bath at a given time, it cannot 
be completely decisive for an actual experiment. The reason is that the rate of change of 
the environment states is order of magnitudes faster than the inverse time scales of the 
experiment. Since these time scales in reality are so different, even the notion of initial 
conditions is not well defined. Relevant initial conditions can be meaningfully prescribed, 
for two interacting systems, only when their intrinsic time scales are not orders of magnitudes 
apart. A theory, which is being introduced to study a real experiment, has to incorporate 
these facts into its structure. It is manifested here by assuming that the reservoir is in 
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thermal equilibrium. This has nothing to do with the consistency of quantum theory, but 
rather with the role of the real initial conditions of the entity we set to investigate. It is 
our opinion that a theory, being a logical construction in terms of a mathematical scheme, 
is supposed to be correlated with the intrinsic approximate nature of the results of real 
experiments, and it has to reflect this facet of physical reality 

The outline of the paper is as follows: In Section II we introduce the model to be in- 
vestigated in terms of the Hamiltonians for the spin-half atom, the massive pointer, the gas 
reservoir, and the interactions among them. Section III is concerned with the motion of the 
atom and the pointer without the reservoir. The quantum interference between the two final 
positions of the pointer is studied in detail. In Section IV we add the bath into the develop- 
ment in time, introduce the Fokker-Planck equation for the pointer in the gas environment, 
and show how the decoherence is set in. In Section V, we discuss the conclusions. The 
details of the solution of the Fokker-Planck equation for the Wigner distribution is outlined 
in Appendix A. In Appendix B we study an intuitive model for the bath, which is simulated 
by a random field perturbing the meter, and compare it with the results obtained for the 
gas reservoir. 



We consider the measurement of the spin components of a spin half atom. We introduce, 
for our microscopic system S, i.e., the spin half atom, the Hamiltonian, Hs, which, in the 
absence of an external magnetic field, has degenerate energy eigenstates 



These are also eigenstates of a^, the component of the spin operator along the 2;-axis, i.e., 
\ cr) — a \ a), where cr = ±1, for the up or down states. We shall take Eq — 0. 
Our measuring device, the pointer, is taken to be a massive particle of mass M (possibly 

a heavy atom), whose center of mass position and momentum are R and P respectively. 

The Hamiltonian, Hp, of the pointer is then 



II. MODEL HAMILTONIAN 




(1) 
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Hp = . (2) 

The pointer is the meter which is assigned to perform a measurement of the z— components of 
the spin of the atom. A measurement is thought of as a process, which generates a correlation 
between the measured property of the microscopic system, i.e. the spin component, and say, 
the final macroscopic position of the pointer. To measure the spin variable a^, which has 



the observable outcomes ±1, we adopt, after Peres [|15|, the following coupling Hamiltonian: 



Hsp = V{t)P,a,, (3) 

where V{t) is a c-number time dependent "impact" function, which couples the spin operator 
to the pointer. The function V{t) has the dimension of velocity, and is active for a very 
short finite interaction time, < t < T. We wish to point out that one massive atom can 
play the role of both the pointer and the spin half system, when the internal microscopic 
variable - the spin - is being coupled for a short time with the center of mass of the atom, 
as in a Stern-Gerlach device for example. 

The effect of the environment on the pointer is to be studied here assuming the following 
picture. The pointer is considered to be immersed in a bath, which is viewed as a set 
of independent identical "field" particles, whose i-th particle is described by the mass m, 
coordinate rj, and momentum pj, and the Hamiltonian 

H^ = y^. (4) 

i 

Each of these bath particles interacts with the pointer via a potential 0(R — rj), and thus 
the pointer - bath interaction Hamiltonian is 

HPB= 5^0(R-r,). (5) 

i 

The total Hamiltonian for our system is thus 

H = Hs + Hp + Hsp + Hb + HpB. (6) 
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The outline of the procedure employed for the measurement of the spin of the atom is 
as follows. Suppose that the atom is initially in the state 

1^5(0)) = 5^a^|a) = a+i |+1) + a_i |-1) , (7) 

where a^- is the amplitude to find the spin in the a state, a* is the complex conjugate of 
Qa, and 0^0(7 = 1- This state is to be measured by the position of the pointer after it 
interacts with the atom. The initial state of the pointer, which is also described by quantum 
mechanics, is assumed to be a Gaussian wave packet: 

V'(R,0) = (27rA2)-=^/4e-^'/^^', (8) 

where A is the spread of the wave packet in R, namely, the initial expectation value of R is 
0, and that of R^ is A^. We take A to be as small as possible, to have a better resolution 
of the position of the pointer. The change of the state of the pointer in time, due to its 
interaction with the atom and the environment, is now studied in details, by solving the 
equation of motion of the pointer and eliminating the degrees of freedom of the "gas" of 
field particles surrounding the pointer. But before we will do that, let us first take a look at 
the case where there is no bath and the measurement of the spin is attempted with a " free" 
pointer. 

III. MEASURING THE SPIN BY A FREE POINTER 

First we consider our system without the environment, and find the time development 
of the pointer from its given initial conditions. Since the "free" spin term of the atom is 
irrelevant, the Hamiltonian is 

+ V{t)P.a,. (9) 

The state of the combined system, atom and pointer, can be given at the time t, by the 
wave function ^(R, (T, t), which is the solution of the Schrodinger's equation 
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with the initial condition 

^(R,a,0) = a. (27rA2)-3/4 e-i?V4A=^. (11) 

It is convenient to introduce the momentum representation for the pointer, namely P |k) = 
/ik |k) , and apply the Fourier transform to the wave function, 



-oo 



^(R,(T,t) = / rf^fc e+*'^-^^(k,(T,t). (12) 



oo 



(27r 

The equation of motion takes the form 

z^-v[/(k, a, t) = (-^ + Vmhaj vl/(k, a, t), (13) 

where k = |k| , and the initial condition is given by 

^(k, a, 0) = a„{8nAY^ e-^'^\ (14) 

The solution of Eq. (|^) is then 

^(k, a, t) = ^(k, a, o)e-^^**-^^W (15) 

where ujk = |g, and X{t) = /J dt'V{t'). 

It is convenient to discuss the measurements in terms of the density matrix of the com- 
bined system. In the (k, a) representation we write the density matrix as 

p(k, a; k', a'; t) = vl/*(k', a', t) v^(k, a, t), (16) 

where ^* is the complex conjugate of From Eq. ( [I5l) we get 



p(k, a; k', a'; t) = p(k, a; k', a'; 0) 

Xg-*('^fc-'^fe')tg-*^W (ka^er-k'^a') ^ (-|^j) 



where, using Eq. (14), the density matrix is initially 
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p(k, a- k', a'; 0) = <,a. {StiA^/^ ^-A^Cfc^+fc-), (ig) 

We now return to real space, and transform the density matrix to the (R, a) representation 
to find: 

p(R, a; R', a'; t) = al,a,{2TiA''i{t))-^'^ 

xexp (-^^^ {c{m-^x{t)af+am-mtyf)) , (19) 

where C(t) = l + it/rf, ^{t) = C*(^)C(^) = 1 + {t/vf- Here tj = 2MA^/h is a characteristic 
time of free quantum diffusion of the pointer, and x is a unit vector along the x— axis. 
Eq. (|1^) is our result for the density matrix of the "free" pointer. In passing we notice that 
we can express formally the density matrix at time t in terms of the density matrix at t = 0, 
namely 

/ + 00 f + OO 

d^Y / d^Y' JiR,R',t;Y,Y',0) p(Y,a;Y',a';0), (20) 
-oo J —oo 

where the propagation kernel is 

J(R, R', t; Y, Y', 0) = / d'k' e'^'^' ^^'-^') / d'k e'M^-^)^~^{^.-^y)t-^x(t) (k.a- 

\ 3 

^ > 'giJ\/[R-Y-xX(t)(T]^/2mg-jAf [R'-Y'-xX(t)cr']^/2m 



We also wish to point out that we could have obtained Eq. ([T9| ) directly from Eq. (|15| 
using 



1 /■+00 

vl/(R,a,t) = ^'fc e-^'^'^'v^(k,a,t) 

= a^{2^A\\t)r^l' 



("4A^f^ ~ ^^W^]') ' (22) 

and p(R, a; R', a'\ t) = ^*(R', a', t)^(R, a, t). 

We now turn to study Eq. ([T9|) . First we observe that our density matrix represents a 
j^nre siate. It is just the evolution in time, by a unitary transformation of the initial product 
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state, ^(0) = (a+i + a_i)'0(R, 0), into an entangled state of the particle and the pointer, 
namely 



are two displaced wave packets of the pointer. In this way the amplitude of the position 
of the pointer is correlated to the spin of the particle. This is not yet a measurement, but 
rather an establishment of quantum mechanical correlations, due to the interaction between 
the "free" pointer and the particle. 

The outcome of a measurement of the position of the pointer is then predicted, using 
the combined system's density matrix, by calculating the probability density to detect the 
center of mass of the pointer around R, i.e., P(R, t) = p(R, cr; R, cr'; i). Since the 

measurement is done after the interaction was completed, namely at a time i > T, we set 
X{t) — X — Jq dt'V{t'), and write P(R, t) as a sum of three terms: 



= a+i^+(R,t) + a_i^_(R, t), 



(23) 



where 



V'±(R) = (2nA\\t))-'/' 




(24) 



P{R,t) = P+i (R, t) + P_i (R, t) + P+1,-1 (R, t) . 



(25) 



Using Aj{t) = A2 



i{t), the first term, with cr = u' = +1, is 



P+i(R,t) = |a+i 1^ {2nA}{t))-'/' 




(26) 



and second term, with a — a' — —1, is 



P_i(R,i) = I a_i 



2 (27rAj(i))-=^/2 




(27) 
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corresponding to the probabilities of up and down spins. While the third term, with a±i 



|a±i|e*'^±, is 



P+i,_i(R,t) = 2v/P+(R,t)P_(R,t) 

XX t 



X cos I --^277\~ + ^~ - ^+ \ ^ (28) 



is due to quantum interference between the two spin states. In Eqs. (|2B|), ( P?] ) and (^) 

Aj(t)=A2[l + (t/r;)2] (29) 

reflects the spatial broadening of the pointer's position as a result of the free quantum 
diffusion. At this point we note that if the "pointer" is a Silver atom, with a mass of 
M = 1.8 X 10^^^ g, and initial spread of position A = 1 fim, 

2MA2 ^ ^ 

Tf = ~ 3 ms (30) 

which is a long time on the scale of a ballistic experiment of an atom. 

The probability distribution for the position of the pointer, long after it departed from 
the "measured" particle, has three components: two of them are "normal", and the third 
one is "strange". As it is clear from Eqs. (p6D and (0), one component represents a positive 
deflection of the pointer by +X, and displays a Gaussian centered around {+X, 0, 0), which 
is correlated with the spin up state. The other one indicates a negative deflection by —X, 
and displays a Gaussian centered around {—X, 0, 0), which is correlated with the spin down 
state of the measured system. The width of each term is of the order of the spread of the 
initial distribution. However the weight of each peek, i.e., the probability of positive or 
negative deflection, is decided by the initial probability of the two spin states, namely by 
I p. The probability distribution of Eq. (|28|) exhibits an interference between the two 
displaced Gaussians. Notice that this interference term oscillates in time with a frequency 
given by 

_ XX 1 
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This interference is an unambiguous outcome of the quantum mechanical development in 
time of the pure entangled state of the pointer and the system to be measured. This result 
is completely reversible in time, and in principle can be reversed to the initial setup. 



IV. THE POINTER IN A BATH 



We turn now to consider the case in which, in addition to the "external" interaction 
with the spin, the pointer is coupled to a reservoir. Our task is to describe the evolution in 
time of the pointer's state under the influence of this reservoir. To this end, we will drive 
a master equation for the reduced density matrix p{t) of the pointer, following the method 



pointer on the bath is assumed to be very small, and the interaction Hps is considered as a 
weak perturbation. The equation of motion for the density matrix of the combined system, 
pointer and bath, is written up to second order in HpB- Finally, by tracing over the degrees 
of freedom of the bath, one obtains an equation of motion for the pointer only, which is 
valid on time scales much larger compared with the typical correlation times of the bath. 
We choose here a particular model for the bath, which is constituted of a set of independent 
field particles, and trace over the dynamical variables of these particles. This results in a 
Fokker - Planck equation for the reduced density matrix, which can then be solved for the 
time dependent behavior of the pointer. 



In order to study the effect of the interaction. Hps, on the pointer, it will be convenient 
to introduce the Fourier transform of the interaction potential, 0(r), in space, in a volume 



outlined in Ref. |T^. The main ingredients of this method are as follows: The effect of the 



A. The Reduced Density Matrix 



V: 




(32) 



q 



Then the interaction Hamiltonian can be written as 
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q 

where the Fourier transform of the density (operator) of the field particles is 

r^(q) = 5^e-'^■^^ (34) 

i 

To get the equation of motion for the reduced density matrix for the pointer, we transform 
to the Interaction Representation (IR), and write the Hamiltonian as 



q 

where, for the pointer 
and for the bath operator 

n(q,t) = e^^^*/^n(q) e"^^^*/^. (37) 

The reduced density matrix of the pointer, in the IR, 

p(t) = TrB {e'^^p+^BWh p(t)e-^(^^p+^s)t/?ij ^ (38) 



obeys the following equation, |T^ Eq. (B.30): 



^ - dti dt2-J2 E <^(q ) 

Trb{[ n(q,ti) A^{ti), [ n(q',t2) Aq'(i^2), W) ® P^]}, (39) 



where the trace is taken over the field particles' degrees of freedom, and is the reduced 



density matrix of the bath. This master equation, Eq. (|39|) is similar to Eq.(l) of Paz and 
Zurek and to that of Unruh and Zurek [0; however unlike the scalar potential that 
simulates their bath, here the environment is a gas. Another general master equation was 



derived by Joos ]T8[ following Pauli, without specifying the bath. 
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Assuming that the distribution of the field particles is stationary in time, and homoge- 
neous in space, we can express the trace of the RHS of Eg. (p9D in terms of the density - 
density correlation function ( CF) of the bath, 

TrB{ Pb n{q,t + r) n(q',t)} = ( n{q,t + r) n{c^,t)) 

(7q(r) = (n(q,t + r)n(-q,t)). (40) 



Eq. (H) using (|g) is written as 0] Eq. (B.33). 

Ap(t) 1 1 V- m2 , 1 f'^^' ,/ 

x{g^{T){A^{t') A^^it' - r) p{t) - A_^{t' - r) p{t) A^{t')) 

+9-ni-r)i m A^^it' - r) A^{t') - A^{t') p{t) A_^{t' - r))}. (41) 

We further assume that the reservoir is in thermal equilibrium in temperature T, and 
thus Pb = e~^^^ /Z, where /3 = 1/kBT is the inverse temperature, fc^ is the Boltzmann 
factor, and Z is the canonical partition function. Introducing the Fourier transform in time 
of the CF, 

1 r"^ 

we can write the spectral density of the CF as 

^?qM = ^5^e-'^^''5^(/i|n(q)|z/)(z/|n(-q)|^) 

x6{E^- + hu), (43) 

in terms of the eigenstates of Hb, i.e., Hb \p) = \p). First we Notice that since n(— q) 
is the Hermitian conjugate of r;,(q), g^iuj) is a real function of q, tu. Then we observe that 
g^f^{—LJ) = e~^^ g^^iuj), which is a statement of the detail balance induced by the thermal 
reservoir on the states of the pointer. 

We consider a simple model for the bath, namely a gas of non-interacting particles 
in temperature T, and average density no. The mass m of the bath particles is taken to 
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be much smaller than that of the pointer, i.e., m/M <^ 1. To calculate the correlation 
function we first express the density operator, n(q) = X]p'^p-q '^P' terms and Op, 
the creation and annihilation operators of a "field" particle in the momentum eigenstate 
I p), of energy ep = ^^p^/2m, and momentum fvp- Since in the interaction representation 
ap(t) = Op exp (— icUpt), where Cp = /icUp, and in thermal equilibrium (Op.q Op Qp'+q Op') = 
5p/ p_q /p-q(l ± /p), where ± is for either Bose or Fermi particles, and /p stands for Bose- 
Einstein of Fermi-Dirac distributions, the correlation function is then 



Ur) = J2 /p-q(l ± /p) e-(— (44) 
p 

and its spectral density is 

g^{uj) = 271 /p(l ± /p+q) 5((^p+q - Up) - uj) , (45) 

p 

We now return to Eq.(O), and express the reduced density matrix for the pointer in its 
momentum representation. We denote by | k) , the eigenstate of momentum ^k, and energy 
= h'^k'^/2M, and find for the matrix elements of the pointer's operator Aq(t) in the IR 

(k I Aq(t) I k') = (k I e'^^*/^ e^'i-^ e-'^Bt/h I ]^/^ 

= e'""^--'' 5k.,k-q, (46) 

where fik,k' = {E]s_—Ei^/)/h. In the Secular Approximation the equation of motion for the 
reduced density matrix, in the Schrddinger Picture can be written as 

|(k I Pit) I k') - l(k I [H„p{t)] I k') = Jk,k'(t). (47) 
Here the effect of the environment is given by the collision term 

q 

X{(G'q(a,k-q) + G';(a',k'-q))(k | p{t) \ k') 

-(G'q(l^k+q,k) + G';(^]k'+q,k'))(k + q | p{t) \ k' + q) } , (48) 

where 
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G^iuj) = / dre'-^ g^ir), (49) 





is the one-sided Fourier transform of the correlation function. If we discard the frequency 
shifts due to the bath, we can replace in the above equation G'q(co') by the real function 
5fq(co')/2. To further facilitate the p— summation in the gq{uj)'s, we make use of the 6— 
functions, and write, for example 

p 

x5{p.q-g+(k,q)}, (50) 



where q = q/g, and 



Q,(k,q) = -q.k-|(l±-). (51) 



B. Classical Gas Bath 

We now consider the bath as a classical gas, where the thermal distribution of the 
particles is Maxwellian, namely 

f, = noi-f' e-^\ (52) 

TT 

with a = /2mkBT . In this case, we have /p ^ 1, and if we replace the summation in 
Eq.(^Ol) by integration over p, i.e. write ^ ^ / d^p/{27i)^, we find 

^ ^r,f-'-^ e- (53) 

and similarly for the other Gq 's of Eq. (|^) . The collision term can be then cast into 

,27r 1 [am 



1 v-^ I ,/ \ i22vr 1 [am 



q 

l2 



vr qh? 

x{( e-^^+e^-'i) + e-"<3+('^'''i)) X (k I p{t) I k') 

g-aQi{k,q) ^ g-aQ2_{k',q)^ X (k + q | p{t) \ k' + q)}. (54) 
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This is still an exact result for a classical gas interacting with the pointer. 

We now make use of the two small parameters in the problem, namely, that the mass 
ratio is small, and that the momentum transfer is small, 

1] = m/M < 1 and g < A;. (55) 

First we expand in small t]. From Eq. (^) we have Q^(k, q) =^ ^ ~''?(k-q=F^) + 0(77^), 
and thus, to first order in rj we can write 

g-aQi(k,q) ^g-a^(l + ^^(k.q^^))_ (56) 

Then we expand (k + q| p{t) |k' + q) in power series of q, to obtain, up to second order in 
q/k 

(k + q| Pit) |k' + q) = (k| Pit) |k') + q ■ + ^) (k| Pit) |k') 

>From the structure of this expansion we realize that it is convenient to express the den- 
sity matrix in terms of the vectors K = (k + k')/2, and p = (k — k') , namely write 
(k| p(t) |k') = p(k, k',t) = p(K,p, t). Substituting these expressions in Eq.(^), assuming 
that 0(q) is isotropic in q, and integrating over q, the collision term is simply 

4,k'W=7^-(K p(K,p,()) 

Here 

plays the role of inverse relaxation time, and D = 'j/2ari plays the role of a diffusion 
coefficient in k space. Using [|19| we can relate the A;— diffusion coefficient, D, in Eq. ( ^81) to 
the standard coefficient of spatial diffusion, Dc = ksT/M'j, i.e., we have D = iM'j/h)'^Dc. 
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Notice that it is spatial diffusion coefficient, Dc, which is connected by the Einstein relation 
to 7. 

The equation of motion of density matrix of the pointer, Eq. (^Tj), becomes, with Eq. (p8[), 
a Fokker - Planck (FP) equation for (k | p(t) \ k'), namely 

I (k| p{t) |k') - i (k| [Hp, p{t)] |k') = 7 • (K (k| p{t) |k')) 

This equation is much like the FP equation derived by Caldeira and Lggett for a harmonic 
oscillator interacting with a bath of harmonic oscillators in equilibrium. The main difference 
is that the friction parameter 7, here is given in terms of the properties of the gas surrounding 
the pointer. 

Before we turn to the solution of this equation in time, we make some estimates of the 
Fokker-Planck coefficient, 7. If the interaction potential 0(r) is specified, we can calculate 
explicitly 7 of Eq.(^). Let us propose a particular interaction potential that will simulate 
the collisions between the pointer and the gas particles. We take 

0(r) = 0oe-^^'/'^^ (61) 

where a is the length scale of the pointer, and 00 is the potential strength. Using this 
potential, the term 10(^)1^ in Eq. (^) can be specified as 



TT^aVoe-"'^'/'. (62) 



For the inverse relaxation time we get 

where g = 2a} j a. Recalling that a = h"^ /2mkBT, and noting that even for objects with 
atomic dimensions ^) ^ 1, we can write 7 as 

7 = ^/5w%(^) , (64) 
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where er = 3fcsT/2 is the average thermal energy and v = ^J^ik^Tjm is the average velocity 
of the gas particles. Although we have an explicit expression for 7, it is still necessary to 
specify 0o- For air in room temperature, a pointer which is a silver atom, we take a = 
1.75 X 10^^ cm to be the radius of the atom, and M = 1.8 x 10^^^ g. Assuming that the 
barrier 0o is, say 50 times the thermal energy e^, we get a typical 7 ~ 2.5 x lO^s"^. 

C. Solution of the Fokker-Planck Equation for the Pointer in a Bath. 

In order to solve Eq. (^Op for the pointer model outlined in section II, we further make 
two simplifying assumptions. First, we observe that since the deflection of the pointer as a 
result of the interaction with the spin is only along the x-axis, the coordinates y and z are 
of little consequence for our purpose. Hence we will consider only the one dimensional FP 
equation along the x- coordinate, and use, for convenience, x instead of the X— coordinate 
of the pointer. If we set h k, k.' ^ k', K = {k + k')/2, and p = {k — k') we can restrict 
ourselves only to 

{k\p{t)\k')=p{k,k\t)=p{K,p,t). (65) 

Our second assumption is that the interaction of the pointer with the atom occurs in a very 
short time, or takes place before the interaction with the bath is switched on. In other 
words, we start from a density matrix which describes a free pointer after its interaction 
with the spin has already been completed, that is we use Eq. (|T7|) as the initial condition for 
the FP equation for p{K,p,t). It is clear that both assumptions are not necessary in order 
to carry out the calculation. The first one comes to simplify the mathematics by removing 
the irrelevant degrees of freedom, while the second one saves us the need to define V{t) 
explicitly in Eq. (j^). We bare in mind of course, that what we are really interested in here, 
is the decoherence process itself, and not so much the creation of the initial coherent state 
of the pointer. 

Suppose therefore that we know the initial density matrix, p{p,K,0), and we seek the 
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solution of p(p, t) at the time t. The free streaming term of Eq. (|60|) for the pointer in its 
momentum representation is then 

- 4 {k\ [Hp,p{t)] \k') = -\{Ek-Ek,)p{k,k',t) 
in in 

Kv 

= thjj^ p{p,K,t), (66) 



Inserting this back into Eq. ( pOD gives us 



d Kv d 

-p{p, K, t) + ihj^ p{p, K,t) = j-^ {K p{p, K, t)) 



+D^pip,K,t). (67) 



where the initial condition can be written with the aid of (|T^ and ([T8|) as 
p{p, K; a, a', 0) = a:,a,2(27rA2) V^e-^A^^^-^^^C--') 

Since the measurement is done after the interaction was completed, namely at a time t > T, 



we 



set X{t) = X = dt'V{t'). 



The method of solving Eq. ( p7| ) with the initial condition ( |58| ) is outlined in appendix 
B. Here we just quote the result, which is 



p{p,K-a,a'-t) = a:,a,e^' /i![!^g-AV/2-.pX(.+.')/2-Dep^g-.^AXp-«e^^*i.^ 

V + u 

where we have defined 

u{t) = l/ADr]{t), 

C\t) 2 (e^* - 1) 



ri{t) 7(e^* + l)^ 



m 

e(t)=^^J [t-A(t)]. (70) 

Once we have this solution we return to the outcome of a measurement of the position of 
the pointer. In the present one-dimensional case the probability density to detect the center 
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of mass of the pointer around x, is determined, as in the previous "free" pointer case by 
P{x, t) = J2a(7' P(-^^ ^' ^' "^'i namely by the diagonal elements of p in real space. We thus 
have to transform Eq. (^) back to {x, x') space and then to set x = x' in the resulting 
density matrix. This can be done conveniently by noting that the Fourier transform from 
{k, k') space to (x, x') space i.e. 

p(x, x', t) = j ^e-^''^' I ^e^'^piK k', t), (71) 

can be written with new variables, 6 = {x — x') and X = (x + x')/2, and the use of the 
definitions K = {k + k')/2 and p = {k — k'), as 

p{X,9,t) = I ge-^l ^e^^'p{K,p,t). (72) 

Now the probability density will be given simply by P(x, t) = J^aa' Pi-^ = x,a]9 = 0, a'; t). 
We write P{x, t) as a sum of three terms: 

P(x,t) = P+i(x,t) + P_i(x,t) + P+i,_i(x,t). (73) 

where P±i(x,t) = p(x, x, ±1, ±1; t) stand for the diagonal elements, and P+i_i(x,t) = 
p(x, X, +1, —1; t) + p(x, X, —1, +1; t) stands for the off-diagonal elements. 

The first and second terms of Eq. ( [73D are the probabilities of the spin-up and spin-down 
components of the measured atom. Performing the transformation in Eq. (^) we get for 
a' = (T = +1 



P^,(x,t) = K,|%/^^ exp <j ^ (74) 



(x - Xf 



and for a' = a = —1 



The spatial broadening of the pointer's position, which is due to both the free quantum 
diffusion and the effect of the bath, is expressed as 

Alit) = A'Ht) + K{t)], (76) 
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where 

2 



and 



x(t) = l+(-i-) (l-e-^T, (77) 



1 \^4A2D 

nit) ={—) ——filt) (78) 



is expressed in terms of the function 



fix) = 2{x- 2^^-^^ + -^V^^. (79) 

^ ^ * (1 + e-^) J (1 + e-^) ^ ^ 



We now compare these expressions with those obtained previously, Eqs. (p6D, (|27| ) and (p9[), 
in the absence the bath. We observe that the previous broadening of the pointer's position, 
A fit), is replaced by A^(t), where x(t) takes care of both the free quantum diffusion and 
the friction caused by the environment, and ^(t) corresponds to the diffusion afflicted on 
the pointer by the gas surrounding it. The effect of the environment on the probabilities 
P±i(x, t) is just to enlarge the uncertainty of the pointer's position. 



Now let us look at the off-diagonal term of Eq. ([73D. Taking a±i = \a±i \ e**^* we get 



P+i,_i(x, t) = 2 v'P+i(x, t)P-i(x, t) exp |-^f|y 



where K(t) is given by Eq.(|78|). This is the main result of the present paper. The function 
P+i._i(2;, t) represents the probability of finding the system in a state of interference between 
the two possible "classical" outcomes of the measurement. If we compare this result with the 
interference term, Eq. (p8|), of the "free" pointer, we observe two changes: (i) the argument 
of the cos in Eq.(^) is replaced by a more elaborate time dependent expression in Eq. (pO|) , 
yet in the limit 7 — > it reduces to the expression in Eq. ([28|) , and (ii) a new time dependent 
exponential factor appears in Eq.(^). This exponential factor is the important effect of the 
environment on the pointer, and we examine it now more closely. Let us write it as e~^^*\ 
where git) is the decoherence function, and it is given by 
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K{t) 



(81) 



2A2(t) 2A2 [/t(t) + K{t)_ ■ 
Notice that it is not the friction, which is due to energy exchange with the bath, that 
is responsible for the decoherence, but rather the diffusion in momentum space, which is 
originating it. 

For times much smaller than the relaxation time, i.e. 7t << 1, we have 

.W^Z,f^V^|(,<f^2Z,fa^Vt (82) 
and since x(t) is of order 1, 



Remembering that D = M'-fkBT/h'^, we write g{t) as (F't)^, where 



r' 



-2 \ 1/3 



IX^-fkeT ( \ 
3A2 VMA2 ) 



1/3 



(84) 



A4 3M" 

plays the role of a kind of inverse decay time of the interference. It is instructive to compare 
r' of Eq. with the analogous expression of Ref. Eq.(3.5). Their result was obtained 
for an oscillator, with a typical frequency of u^,, which is coupled to a reservoir of harmonic 
oscillators, in the high temperature limit, for the weakly damped case, 'j/ujr ^ 1, and 
at initial times. The only difference between their and that of our system, F', is that 
their oscillator's quanta huj^i is replaced by the energy term h'^/MA'^, which is related to 
the quantum mechanical energy associated with the initial spread of the pointer. To make 
a simple estimate for T' let us take our pointer as a silver atom. We thus take A ~ 1 
fim, to represent a pinhole for localizing the atom, X ~ 1 cm, as a deflection position, 
M ~ 1.8 X 10^^^ g, and a = 1.75 x 10^^ cm. For air at room temperature, namely with 
Tr = 1/7 — 0.4 ns, we obtain F't,. ~ 50. We see that the typical decoherence time 1/F' is 
almost two orders of magnitude smaller then relaxation time. However this behavior does 
not proceed for times much longer than the relaxation time. When time increases, so does 
g{t), and eventually it reaches a saturation value of g{t) ^ ^ (X/A)^ , which is of the order 
of 10^. Nevertheless, at later times when jt > 1, but long before saturation, we have 
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n{t) 2D ^— j t, (85) 



which gives 

2 



9it) 



2A2 



2 



(86) 



2d{^) t+1 

At this time range, g{t) is approximately hnear in time, i.e., g{t) Ft, with a slope given 
by 



X knT X kuT 



(87) 



A4 M7 A2 ;i2^ MA2 ' 
for the silver atom it yields r/7 = 3.5 x 10^. Thus the typical decoherence time l/F in the 
linear part is much shorter then l/F'. The meaning of this result, as was pointed out by Ref. 
0, is that as far as the decoherence is concerned, the system becomes aware of the bath, in 
a time scale much shorter than Xf,, and the interference is deteriorating almost immediately. 
Our decoherence rate, F, of Eq. (|87|), for the gas bath, is quite similar to that of the 
oscillators' bath of Ref. Eq.(3.7), at the high temperature and strongly damped limit. It 
is interesting to note that the strongly damped limit, which is defined as Rose = i/^r ^ 1 
for the oscillators' bath, is replaced here by Rf = 'jTf ^1. In our analysis the meaning 
of the "strongly damped" case is not clear, yet taking h^/MA"^ instead of the oscillator's 
quanta hwR of Ref. 0, we see that Rf = 'yTf ~ 10^ is indeed very large. So in that sense 
we are in the strongly damped limit. It is also interesting to compare F of Eq. (|87|) with 
the decoherence rate, Tz — 7(X/Ar)^, quoted by Zurek Ref. where At = h/ ^/2Mk^r' 
is the thermal de Broglie wave length of the pointer. Apparently Zurek assumes that the 
decoherence time, is proportional to the inverse of the coefficient of the last term in his 
equation (9). We find that the ratio T /V z is of order 1/-Rj, i.e., the decoherence time due 
to the interaction of the pointer with the gas is much longer than Zurek estimated. This 
is also the case for the Caldeira and Leggett P decoherence rates, Fi and F3: The ratio 
Ti/T z in the weakly damped limit, is of order Rocs/ Rf-, while in the strongly damped limit, 
Fa/Fz is of order l/RfRosc- 
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Fig. 1 shows a graphical description summarizing the discussion above. The decoherence 



function of 7t. In the inset, we see more closely the behavior at initial times, where g{t) 
starts as (F't)^ and quickly goes into the linear behavior as Tt. 



In this paper we have studied the affect of the environment on the process of measurement 
of the state of a microscopic spin particle by a classical pointer. It has been shown that, 
the entangled state, which occurred due to the interaction of the spin particle with the 
measuring apparatus, the pointer, is developed into a mixed state by the environment, in an 
extremely short time. An initially prepared state of a spin half particle, as a superposition of 
up and down states, is transformed by the interaction into corresponding two distinguishable 
positions of the pointer. The probability to find the pointer in these positions is expressed in 
terms of two separated wave-packets, and an interference term, which is a manifestation of 
the entanglement between the particle and the pointer. Due to the coupling of the pointer 
with the environment, which is modeled here as a gas of independent particles, rather than 
an ensemble of harmonic oscillators, this quantum interference decays extremely fast in 
time. The final mixed state, which is the outcome of the decoherence, introduced by the 
environment, is being expressed as the two separated wave-packets of the pointer's position. 
The probabilistic nature of quantum mechanics is still maintained in the final representation 
of the pointer, but the interference is destroyed. An exponential decay function, which 
portrays the decoherence, is expressed in terms of the parameters of the pointer and the 
bath. 




normalized to the saturation value 




V. CONCLUSIONS 
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A. Solving the Fokker-Planck Equation for the Pointer 



In the present Appendix we seek the solution for Eq. (|67[) with the initial condition (|68D. 
It will be convenient to study the density matrix as a Wigner distribution function, which 
is defined by the transformation 

/oo 
dX e-'^P W{X,K,t). (88) 

Eq. ( |67D is transformed to the equation 



^W{X,K,t) + ^^W{X,K,t) 



= {KW{X,K,t)) + D^W{X,K,t), (89) 

where K, t) plays a role analogous to that of a classical distribution function, f{x, p, t) 

to find a massive test particle near the position x and the momentum p, while moving in 
a gaseous bath. This partial differential equation can be solved by the use of the method of 
characteristics. We start by ignoring the term with the second order derivative in K, and 
we look for curves that satisfy the conditions 

dX , dK dW 

dt = = -— . (90) 



hK/M -fK -fW 

Integrating with respect to t we have 



dK 

dt = -- = Ke-^' ^ K = Ke^* (91) 

'-)K 



dt-X = -47'^e--* + e-e = X + |^ (92) 



hK/M 7M ^ ^ M7 

dt = ^^W = ue-^' (93) 

where we wrote the constants of integration over t, i.e. k, ^ and uj, in terms of the variables 
K, X, t. 

We can now make the independent variables' transformation K,X,t ^ k, ^, r, where 
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hK 



r = t (94) 
and find that in terms of tlie new variables, and tlie substitution 

Ty(e,K,r) = ^(e,«:,r)e^^ (95) 

Eq. ( ^9|) takes tlie form 



= + — .«...r). (96) 

This equation is a diffusion equation in ^, k, r space, with a time dependent diffusion co- 
efficient. This can be solved for example by using Fourier transform in space, namely, 
introducing f{p,x, r) by 

-I roc roc 

fip, X, r) = — / e-^«P / dx e'^M^, k, t) (97) 

J -co J —oo 

Eq. (RBp is simply replaced by 



with the formal solution 

^ ' -p ) !> , (99) 



/(p, X, r) = /(p, X, 0) exp | ^ (^e^'^'x 
where f{p, x, 0) is the initial value of /. If we denote by 



M7 



r]{t)= / dt'e^^'' 



e 



27* _ 1 



27 

C(t) = / dt'e^' = , (100) 

Jo 7 



fip, X, t) = f{p, X, 0) exp <j -D { r]ir)x' - 2C(r) -^xp + ) |^ • (101) 



we can write Eq. (00) as 
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We return now to Eq. ([97D to find 

-| poo poo 

k,t) = / dp e*P« / dx e-'""" 

271" J —oo J —oo 



1 poo poo poo poo 

= — — ^ / dp e'P^ / dx e"'"" / d^' e-'^'P / dn' e'"'" 

/ J —oo J —oo J —oo J —oo 

^-D{,irW-2ar)^.pM^frp-) ^(^/^ ^' ^ Q). (102) 

Note that at t = we have W{^, k, 0) = u;(^, k, 0) and the transformation in Eq. ( p4|) has 
the form 



M7' 

hence we can write Eq. (|102|) in terms of the Wigner distribution 

-1 roo roo 

'^'^ J —CO J —00 

^-D{vir)x^-2ar)^xpH^?TpHia-a')Xo[viT)(cT-a')Xo-2aT)^p+2viT)x]) j^^^^ 

1 /'OO /»oo 

= — ^e^* -.Xe^*(<x-.')Xo / dp ^ipi^+m / rfx 
(27r)^ J_oo J-00 

rfX'e^H^'+wjM dirV^'" iy(X',ir',0). (103) 



00 



Using Eq. (|94D and Eq. (|95D the Wigner distribution at time t can be then expressed in 
terms of its initial value as 

1 poo poo 

W{X,K,t) = — dX' / dK' J{X,K,t;X',K',0) W{X',K',0), (104) 
'^'^ J-00 J —00 

where J is the propagator function for the Wigner distribution, 

-| poo POO 

J{X, K, t- X\ K', 0) = e^* - / dp [(^-^')+if,(i^-i^')] / ^x e-(<^^*^-^') 

2'''^ J -00 J-00 

^_D{vit)x^-2at)j^xpHj^rtp^) _ ^^Q5^ 

We can write this propagator explicitly by preforming the integrations over x and p, which 
gives 
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J(X, K, t- X', K\ 0) = e^* ^j^^-uie-^'K-K'f^-'^/iDe^ ^-^^^^ 



where we used the following definitions 



u{t) = l/ADr]{t) 



riit) 7 (eT* + 1) ' 
^{K,K\t)^^m{K + K'\ 

^X,K,X',K\t) = [{X-X')-'d]\ 



Since we are interested in the density matrix in {K,p), we proceed by transforming back 
the Wigner distribution in Eq. (|104|) using the inverse of Eq. (|88D , and write 

p{p,K;a,a';t) = — dX' dK' W{X' , K'; a, a';0) dX e-'^'' J{X, K,t; X' , K' ,0) 

J —CO J — OO J —CO 

1 /'CXD pCO ^ 

= — / dK' \ dX' J{p,K,t;X',K',0)W{X\K';a,a';0), (108) 

V-OO ^— CJO 

where the transform of the propagator is given explicitly, with the help of Eq. ( |106| ), by 

/OO 
dX e-'^PJ{X, K,t;X',K',0) 
-OO 

= e^* v^e-"(^'*^-^')'e-*^'fe-^®^''-*'^f. (109) 



Looking at the integral over X' in Eq. ( |108| ), we note that the X' dependence of the 
propagator J{p, K, t; X', K', 0) as it appears in Eq. (|109|) , is only through the factor e"*^'^. 



Hence this integral is actually a Furrier transform of the initial Wigner function, which bring 
us back to the initial density matrix 

/OO 
dX' e~'^'m[X\ K'] a, a'; 0). (110) 
-OO 

This last result enable us to write the density matrix at any time t, in terms of the density 
matrix at time t = 

pip,K;a,a';t) = — / dK' J^ip, K,t; K\0)pip, K'; a, a';0), (111) 
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where here the propagator is given simply by Eq. (|109|) without the e factor 

J,(p, K, t; K\ 0) = e^* v^e-"(^'*^-^')'e-^®p'-*'^^ (112) 



Consider now the initial condition for the density matrix (68 



p(p, K- a, a', 0) = a;,a.2(27rA2) V2e-2A^^^-^^o(<x-.')e-AV/2-ipXo(.+<x')/2_ ^^^g^ 

Substituting Eq. ( |113| ) into Eq. (|1 1 1|) , using Eq. ( |112| ) and preforming the integration over 
K', we get the explicit form of the density matrix in {p, K) space, at time t, 



' 2A^ + u 

xexp^-[-^'°'---'>;^^^t""'"''^'' h (114) 

4(2A2 + m) ^ ^ 

which is our goal in this derivation. 



VII. APPENDIX B 

A. The Pointer in a Random Velocity Field 

In Section IV we derived the master equation for the behavior of the pointer under the 
influence of the reservoir. Our starting point was the explicit form of the interaction Hamil- 
tonian 0, which after the averaging over the bath determined the inverse relaxation time 7. 
As a complement to this explicit derivation, it will be instructive to look at a much simpler 
model, in which the effect of the environment is described in terms of a random velocity 
field, v(t). The interaction of the pointer with the bath, which is simulated by this random 
field, perturbs the meter and introduces decoherence. Namely it destroys the interference 
by introducing randomization of the phases, and making the pointer's readings distinguish- 
able. The development in time of the entire system is described by unitary transformation, 
which originates probability into the dynamics by quantum mechanics. This simple model 
is interesting because in addition to being readily understood, it gives essentially the same 
behavior as the more elaborate treatment of section IV. 
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Thus taking the Hamiltonian (|^) and adding the random field interaction we have 

H=^ + Vima. + y^it)-P. (115) 

The effect of the environment on the pointer is assumed to be represented by the interaction 
Hamiltonian 

HpB = v(t) ■ p. (116) 

Here v(t) is a random vector field, which simulates the random impacts suffered by the 
center of mass of the pointer, while in contact with the bath around it. The reservoir 
itself can be thought of as an ensemble of many external particles, which collide with the 
pointer. The choice of this Hamiltonian can be thought of as an extension of the spin particle 
Hamiltonian, Hsp- Namely, it is viewed as corresponding to many individual impacts on the 
pointer momentum P, much like the impact on the center of mass, due to the atomic system, 
Eq. (^). In our simple model the reservoir is not represented by a Hamiltonian, but rather 
by the characterization of the random "velocity" field, to be introduced at the proper stage 
of the analysis. 

This turns Eq. ( |lOl) into 

«*(R.-.<) = (-^^-«v-(tvA 

-iSv(0-^)*(R,ff,i), (117) 

and Eq. (0) into 

^^^(k, a, t) = + hV{t)k.,a + m<. ■ v(t)^ ^(k, a, t). (118) 

The solution of this equation is then 

^(k, a, t) = ^(k, a, o)e-*^'=*-^^W k^a-ik-^t)^ ^^^g^ 

where 

x(t) = [ dt'v{t'), (120) 
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is also a random field with the dimensions of length, and the density matrix is 
p(k, a; k', a'; t) = p{k, a; k', a'; o)e-^(^'=-^fc')* 

Since x(t) is a random field it should be expressed in terms of a distribution, iy[x(t)]. 
The velocity field is a sum of many "impulses" , so does the length random field, x(t), namely 

i i 

x(t) = 5^ X, = ^ ^rm. (122) 

Here Vj is the "strength" of the i-th velocity impulse, f{t — ti) samples the i-th impulse, 
whose duration is Tj, and the summation is carried out over all the impulses which occurred 
during the integration time t in Eq. ( |120| ). Assume that the probability density of the 
individual Xj is Gaussian, namely that 

P(x,) = {2na^r'/' exp {-^^/2a^) . (123) 

The single parameter, (jj, was introduced to represent the width of the Xj distribution, which 
is also assumed to be spherically symmetric. The probability density to find x(t), at the time 
t, to be between x and x + dx is then 

Wi^t)] = n / ^^^^(^^) '^[^W - J2 ^j]- (124) 

i 3 

We now transform the (5-function as 

1 r+oo 

5[x(t) - = 7^ j ^ ^qe-^^-'^^'^-^^ (125) 

and Eq. ( p.24| ) is written as 

X J] f dx,P(x.) e'l-"'. (126) 
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Performing the Xj-integration using Eq. ( |123| ) 



+00 

2„2 , 



c/xiP(xi)e^'»-''' = e-"^« /^ (127) 

the product of the exponentials yields 

1 r+00 

WW{t)\ = —— t;qe-''i"We-(i/2)5'^'"»'. (128) 

Given that the impulses' rate, i.e., the number of impulses per unit time is u, and that the 
average width of the CTj-distribution is a, we get 

Y^a^^uaH, (129) 

i 

and finally, 

+00 

iq-x(t) -(l/2)q2!ya2t 



1 r+oo 



Once the Distribution of x(t) is known, we can find the average of the density matrix in 
k-space. Notice that we could first return to the real space by Eq. (|T^, and then perform 
the averaging. The average over the random field is simply 



rfxiy(x)e-^(^'^'>^ = e-(^'('^-'^')V2). (131) 

J —00 
where 

^2 = ;^^2^_ (^^32) 



Notice that Eq. ( |131|) is similar to the characteristic function of the distribution in x. 



The density matrix of Eq.( |121|) , averaged over the distribution function of Eq. ( |130|) , and 
denoted by p, is 

p(k, a; k', a'; t) = p(k, a; k', a'; o)e-*(^'=-"^')* 

^g-iX (fe,<7-fcia')-/32(k~k')V2_ (133) 
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We now use Eq. (0) and Eq. (0) to return to real space and find 



p(R, (J] R', a'] t) = <,a^(87rA2)3/2 



(2vr) 



X 



+00 r+00 

'3iJ „-jk'-R' 



^3^ gik-R g-A^(k^+k'^)g-»(c^fc-c.,,)t-jR-(k<7-k'a')-/3^(k-k')^/2^ ^^-^3^^| 



where, for convenience of calculation, we have introduced a deflection vector 



R=(X,0,0). (135) 



We notice that the integral in Eq. ( |134|) is a product of three equivalent integrals. 



p(R, a; R', a'; t) = al,ajjyh, (136) 



where, for i = x,y, z, 



1/2 



exp 



2nA}{t) 
Here we have introduced 



dk'. e- 
1 



[am? + c*(t)^' + (f^: - a)'(/5V2A^)) 



(137) 



— Ri — RiCT, 



n[ = R!i- Ria' 



A,W = A'[e(t)+(/97A')] 



(138) 



where in the last equation, A^(t) reflects the spatial broadening of the pointer's position. 
This is due to both the free quantum diffusion, and an additional diffusion term, P"^, caused 
by the random field v(t), which simulates the environment. Using Eq. (|137|) we can write 
the three dimensional density matrix of Eq. (|134| ) as 

3/2 

p(R, a; R', a'; t) = a^itta 



2vrA2(t) 



X exp 



4A2(t) 



(139) 
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where we now have 



12 = R - xX(T, 

n' = R' - xXct'. (140) 



Eq. (|139|) is our result for the reduced density matrix of the pointer in space, after the 



environment was "traced out". Notice that the spin indices can be considered as parameters 
for the pointer's position. It is clearly seen that when (3 —>■ Eq. ( |139| ) is reduced back to 
Eq. (0). 

We now turn to calculate the probability to locate the pointer in the position R in space. 



We set R' = R in Eq.( p,39| ), and find the diagonal elements in space of the density matrix, 
p(R, cr; R, a'; t). Again we write this spatial probability density, as in Eq. (pSj) , as a sum of 
three terms. The first term is correlated with the spin up state, and we set a' = cr = +1, or, 

fi' = fi = R - xX, (141) 

which leads to: 



3/2 



P+i(R,t) = I a+i 



2nA}it) 
1 



^^^P|-^|R-x.Y]^j. (142) 

The second term is correlated with the spin down state, and we set a' = a = —1, or, 

O' = O = R + xX, (143) 



and find 



3/2 



P_i(R,t) = I a_i 



2nAl{t) 



xexp|-^^[R + xX]^l. (144) 

The probabilities of Eqs. ( |142| ) and (|144| ) are similar to these of Eqs. ( p6D and (0). The 
only difference is that the free diffusion spread of the Gaussians, Aj(t), is replaced by A^(t) 
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of Eq. (|138|) , which adds the effect of the environment on the pointer's position. We thus 



observe that the probabihty to find the pointer, either deflected by +X for the spin up case, 
or by —X for the spin down case, is "performing" a diffusion-hke motion in space, which 
is strongly affected by the random impulses inflicted by the environment. It is interesting 
to notice that even in the absence of interaction with the spin particle (i.e. if X = 0), the 
pointer's position would have shown diffusion-like motion of its Gaussian wave-packet, with 
the same time dependent spread, A^(t). 

This diffusion-like behavior of the pointer, due to the interaction with the environment. 



is a reminiscent of the classical diffusion of a Brownian particle immersed in a fluid |]l9l. If 
we compare the result of Eq. (|130|) with the analogous one of the Brownian motion, we find 
that the diffusion coefficient is D = /?^/(2t) or with Eq. ( |132|) D = z/a^/2. This can provide 



us with an estimate for a. Assuming that the environment is a bath in thermal equilibrium 
with temperature T, then, in analogy with Brownian particle, the classical limit yields 
2kBT/{Mu'^), where ks is Boltzmann's constant. If we consider the environment as 
an ensemble of harmonic oscillators, with a typical frequency of ujq, we get 



( 2 \ ^0 



which approaches the classical limit when {huo/kBT) <^ 1. The collision rate, z/, which 
is the inverse relaxation time, Tr, of the pointer, depends of course on the character of the 
bath. 

It is also instructive to estimate the time scale, on which this diffusion broadening be- 
comes so wide, as to erase the record of the measurement. We can estimate this time, tbiuer, 
by noting that when the width Ap{thiuer) is comparable with the deflection of the pointer 
X, it becomes difficult to distinguish between the spin up and spin down results, which 
were previously registered. Thus taking X /A^ ^ 1, which is a necessary condition for a 
meaningful pointer, and ^(t) ~ 1, we have 

tbiuer ^ r,{X/af. (145) 
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This "theoretical" time of bluer is a consequence of the classical Brownian-hke diffusion of 
the pointer, and since always X ^ a, this time is much greater then the relaxation time. 

We turn, now, to the off-diagonal term of the spatial probability of the pointer, namely 
to the interference term, which is a purely quantum mechanical effect, like that of Eq. (p8|). 



We set, again, in Eq. ( |139| ), R' = R, but now take a' = —1, and a = +1, and have 

n' = R + xX, 

rj = R - xX, (146) 



thus using Eqs. (|142| ) and ( |144|) we get an expression similar to that of Eq. (pHD, i.e.. 



P_i,+i(R,t) = 2y/p+{R,t)P4R,t)cos I ^1 + - ) exp ' ^ 



(147) 



Comparing this interference probability with that of the free pointer, Eq. (p8|), we notice 
that in Eq. ( |147| ), beside the replacement of Aj(t) by A^(t), an additional exponential factor, 
e"^^*-*, appears, with 

This damping effect stems, uniquely, from the interaction with the bath, and it is similar 
to that found in Refs. This exponential factor will kill the interference term almost 

instantaneously. To show it, we rewrite the damping function, g{t), of Eq. ( |147D as 

where we took ^(t) = 1. First we investigate g{t), at early times, when t <^ Tr{^y. We see 
that g{t) —>■ t/Tint , where the interference damping time is given by 



Tint =2rr\^ = j , (150) 

which is indeed extremely small compared to r^. To make a simple estimate let us take 
the relevant parameters for a silver atom pointer i.e. A ~ 1 /im, a ~ 0.1 /im and X ~ 1 
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cm.This yields Tint /'Tr — 10"^'^. It is interesting to compare Eq. (|149| ) with Eq. (|86D 
of the gas reservoir modeL We see that both expressions are formally identical. Taking 
D = '-yMhsT 1%^ and tj = 2MA'^/h we can write the term containing t in the nominator of 
Eq. as 

The analogous term in Eq. (|149D gives 



where we have taken a"^ —>■ 2kBT /{Mv'^), and = 7 = l/r^. We notice that in the regime 
where g{t) is linear in time both the Fokker - Planck and the random field models give the 
same behavior. Of course when the time t increases, so does g{t), and eventually it reaches 
a saturation value of g{t) — > \ (X/A) . 
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VIII. FIGURE CAPTIONS 



Fig. 1 

A graphical description of the decoherence function, g(t), in the exponent of Eq. 
normahzed to the saturation value | (X/A) and as a function of •yt. In the inset, we see 
more closely the behavior at initial times, where g{t) starts as (T't)^ and quickly goes to the 
linear behavior as Tt. 
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